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A  two-dimensional  model  of  a  flowing-electrolyte  direct  methanol  fuel  cell  has  been  developed  to  predict 
the  performance  of  the  cell  under  various  operating  conditions.  Governing  equations  including  the  pro¬ 
ton  and  electron  transport,  continuity,  momentum,  species  transport  for  methanol,  water,  and  oxygen, 
and  the  auxiliary  equations  are  coupled  to  determine  the  output  parameters.  These  parameters  are  the 
concentration  distribution  of  the  species,  cell  voltage,  power  density,  and  the  electrical  efficiency  of  the 
cell.  After  validation  with  the  experimental  data,  several  simulations  are  carried  out  to  study  the  effects 
of  the  fluid  velocity  at  the  fuel,  air,  and  flowing  electrolyte  channel  inlets  on  the  output  parameters.  In 
addition,  the  effect  of  recirculating  the  methanol  at  the  flowing  electrolyte  channel  outlet  is  assessed.  The 
results  show  that  higher  fluid  velocities  at  the  fuel,  air,  and  flowing  electrolyte  channel  inlets  are  needed 
to  obtain  higher  power  densities.  However,  an  increase  in  the  fluid  velocity  at  the  fuel  channel  inlet 
causes  a  decrease  in  the  electrical  efficiency  of  the  cell.  It  is  also  found  that  the  electrical  efficiency  of  the 
FE-DMFC  can  be  further  increased  if  the  methanol  leaving  the  flowing  electrolyte  channel  is  recirculated 
into  the  methanol  storage  tank. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Direct  methanol  fuel  cell  (DMFC)  is  one  of  the  most  promising 
technologies  that  can  be  used  in  portable  power  applications.  This 
fuel  cell  is  especially  crucial  for  devices  where  high  power  den¬ 
sity  and  long  operation  time  are  required,  e.g.  laptops  and  mobile 
phones  [1].  The  main  distinctive  feature  of  DMFC  is  its  opera¬ 
tion  with  liquid  methanol.  Some  of  the  advantages  of  the  liquid 
methanol  are:  easy  to  store,  low  cost,  and  high  energy  density. 
However,  the  cost  of  some  of  the  components  such  as  catalysts 
is  still  a  major  issue  for  its  successful  commercialization.  In  addi¬ 
tion  to  its  expensive  components,  main  technical  challenges  with 
DMFC  include  the  low  power  density  and  electrical  efficiency  due 
to  the  methanol  crossover  problem  and  the  low  rate  of  methanol 
oxidation  kinetics  on  the  anode  [1,2]. 

There  has  been  a  significant  progress  in  the  DMFC  research  in 
last  decade  to  overcome  the  technical  issues  mentioned  above  and 
to  increase  its  performance.  Controlling  the  operating  parameters  is 
one  of  the  methods  used  in  this  regard.  For  example,  the  most  well- 
known  approach  to  minimize  the  effects  of  methanol  crossover  is  to 
limit  the  methanol  concentration  at  the  fuel  cartridge  to  levels  such 
as  5  wt%  [3].  Choosing  appropriate  materials  for  the  components  of 
the  cell  also  plays  an  important  role  for  the  improvement  of  the 
performance.  The  common  material  for  the  membrane  is  Nation®. 
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This  material  is  relatively  durable  and  has  high  ionic  conduc¬ 
tivity  and  chemical  stability;  however  it  has  high  methanol 
and  ruthenium  crossover  rate,  high  cost,  low  temperature  limit, 
and  high  humidification  [4].  Neburchilov  et  al.  [4]  reviewed 
alternative  materials  that  can  be  used  instead  of  Nation®  in  a 
DMFC.  They  conclude  that  hydrocarbon  and  composite  fluori- 
nated  membranes  currently  show  the  most  potential  for  low  cost 
membranes  with  low  methanol  permeability  and  high  durabil¬ 
ity.  Platinum-Ruthenium  (Pt-Ru)  is  commonly  used  at  the  anode 
catalyst.  As  the  usage  of  Pt  increases  the  cost  of  the  system  sig¬ 
nificantly,  Serov  and  Kwak  [5]  reviewed  the  non-platinum  anode 
catalysts  for  DMFC  applications.  Their  study  showed  that  the  per¬ 
formance  of  non-platinic  catalysts  is  much  lower  than  that  of 
the  platinum.  However,  new  catalysts  based  on  metal  free  carbon 
nitride  nanotubes  are  promising  support  material  for  more  active 
anode  materials.  Kamarudin  et  al.  [6]  discussed  in  their  review 
paper  that  the  combination  of  DMFC  with  thin  film  batteries  (i.e.  a 
hybrid  power  system)  is  one  of  the  possible  short-term  solutions 
to  overcome  the  economical  issues  associated  with  DMFC.  In  addi¬ 
tion  to  controlling  the  operating  parameters  and  the  search  in  the 
alternative  materials,  new  configurations  and  designs  have  been 
also  proposed  to  increase  the  performance  of  the  DMFC.  Among 
them,  one  of  the  most  promising  designs  is  the  Flowing-Electrolyte 
Direct  Methanol  Fuel  Cell  (FE-DMFC)  [7],  which  is  discussed  below 
in  detail. 

The  FE-DMFC,  which  was  proposed  by  Kordesch  et  al.  [7],  is  a 
novel  DMFC  design  which  provides  performance  improvement  by 
eliminating  the  methanol  crossover  from  the  anode  to  the  cathode. 
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Nomenclature 

ciio 

exchange  current  density  times  specific  area 

(Am-3) 

C 

concentration  (mol  m-3) 

D 

coefficient  of  diffusion  (m2  s-1 ) 

EEq 

equilibrium  voltage  (V) 

F 

Faraday  constant  (sAmol-1) 

i 

current  density  (Am-2) 

j 

volumetric  current  density  (Am-3) 

Jxover 

crossover  current  density  (A  m-3 ) 

Kc 

methanol  oxidation  reaction  constant  (molm-3) 

l 

length  (m) 

LHV 

lower  heating  value  in  molar  basis  (J  mol-1 ) 

MW 

molecular  weight  (g  mol-1 ) 

n 

normal  vector 

nd 

electro-osmotic  drag  coefficient 

O  repeat 

number  of  repeat  elements 

N" 

molar  flow  rate  per  cross  section  (mol  m-2  s-1 ) 

P 

pressure  (N  m-2) 

R 

universal  gas  constant  (J  I<-1  mol-1 ) 

Sc 

source  term  in  the  conservation  of  species 
(molm-3  s-1) 

Su 

source  term  in  the  momentum  equation  (N  m-3) 

t 

thickness  (m) 

T 

temperature  of  the  cell  (K) 

u 

fluid  velocity  (m  s-1 ) 

V 

average  fluid  velocity  (m  s-1 ) 

y cell 

cell  voltage  (V) 

V 

flow  rate  (m3  s-1) 

X 

distance,  m;  molar  ratio 

y 

distance  (m) 

w 

width  (m) 

power  density  of  the  cell  (W  m-2 ) 

Greek  letters 

a 

transfer  coefficient 

£ 

porosity 

P 

density  (gm-3) 

o 

conductivity  (S  m-1 ) 

K 

permeability  (m2) 

1 

stoichiometric  flow  coefficient 

IX 

dynamic  viscosity  (N  m-2  s) 

1 1 

overpotential  (V) 

Pel 

electrical  efficiency  of  the  cell 

0 

potential  (V) 

T 

shear  stress  (N  m-2) 

Subscripts 

a 

anode 

ac 

air  channel 

c 

cathode 

eff 

effective 

fc 

fuel  channel 

fee 

flowing  electrolyte  channel 

in 

inlet 

l 

electrolyte  phase 

MeOH 

methanol 

s 

electrode  phase 

xover 

crossover 

Superscripts 

AC 

air  channel 

FC 

fuel  channel 

in  inlet 

ref  reference 


Schematics  of  a  DMFC  and  a  FE-DMFC  are  shown  in  Figs.  1  and  2, 
respectively.  As  can  be  seen  from  these  figures,  the  FE-DMFC 
requires  additional  layers  such  as  an  additional  membrane  to  sep¬ 
arate  the  flowing  electrolyte  from  the  catalyst  layer  and  a  flowing 
electrolyte  channel  in  which  the  diluted  sulfuric  acid  flows  to  carry 
away  the  crossover  methanol.  This  design  basically  reduces  the 
cathodic  overpotential  as  it  hinders  the  amount  of  methanol  reach¬ 
ing  the  cathode  catalyst  layer.  Flowever,  the  ohmic  losses  increase 
due  to  the  inclusion  of  the  additional  layers.  Diluted  sulfuric  acid  is 
generally  considered  as  the  fluid  to  be  pumped  through  the  flowing 
electrolyte  channel  since  it  has  high  proton  conductivity  to  reduce 
these  additional  ohmic  losses.  In  addition,  the  performance  of  the 
cell  could  be  further  increased  if  the  methanol  recovered  from  the 
outlet  of  the  flowing  electrolyte  channel  could  be  separated  from 
the  sulfuric  acid  and  recirculated  back  to  the  fuel  channel  inlet.  The 
application  areas  of  the  FE-DMFC  are  more  limited  compared  to  the 
DMFC  due  to  the  additional  system  components.  Some  examples 
of  the  areas  that  FE-DMFC  could  be  used  include  backup  power  for 
recreational  activities,  golf  cars,  and  forklifts. 

Performance  of  DMFC  and  FE-DMFC  can  be  estimated  through 
modeling  studies.  There  have  been  numerous  papers  published 
on  DMFC  modeling.  Some  of  the  recent  papers  have  included  the 
multi-dimensional  effects  to  predict  the  distributions  of  the  out¬ 
put  parameters  in  space  (e.g.  [8-11]),  the  two-phase  transport 
taking  into  account  the  effects  of  the  gaseous  C02,  CH3OH,  and 
H20  (e.g.  [10-12]),  and  the  detailed  water  transport  phenomena 
considering  the  combined  effects  of  the  hydraulic  permeation, 
electro-osmotic  drag,  and  diffusion  at  the  thin  membranes  (e.g. 
[13-15]).  As  opposed  to  the  many  efforts  in  DMFC  modeling,  there 
have  been  only  a  few  papers  published  in  FE-DMFC  modeling. 
Kjeang  et  al.  [16,17]  studied  the  methanol  crossover  reduction  in 
the  flowing  electrolyte  channel  with  different  operating  parame¬ 
ters;  however  the  results  of  their  study  were  limited  since  they  did 


Fig.  1.  Schematic  of  a  2D  cross-section  of  the  DMFC. 
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not  include  the  transport  phenomena  in  the  backing  layers  as  well 
as  the  fuel  and  air  channels.  In  a  recent  publication,  Colpan  et  al.  [  1 8  ] 
developed  a  one-dimensional  FE-DMFC  model  to  study  the  effects 
of  flowing  electrolyte  channel  thickness,  concentration  of  methanol 
at  the  inlet,  volumetric  flow  rate  of  the  flowing  electrolyte,  and  the 
recirculation  of  the  methanol  at  the  flowing  electrolyte  channel 
outlet.  Their  study  showed  that  the  maximum  power  densities  of 
the  DMFC  and  FE-DMFC  are  almost  same;  but  the  electrical  effi¬ 
ciency  of  the  recirculated  FE-DMFC  could  be  57%  more  than  that  of 
the  DMFC.  In  another  publication,  Ouellette  et  al.  [19]  developed 
an  one-dimensional  two-phase  model  to  predict  the  effects  of  C02 
formation  on  the  performance  of  the  FE-DMFC.  Their  study  showed 
that  the  presence  of  C02  results  in  a  reduction  in  the  methanol 
crossover  as  well  as  the  active  area  of  the  fuel  cell. 

The  literature  survey  discussed  above  revealed  that  there  has 
not  been  a  paper  published  on  the  multi-dimensional  model  of  a 
FE-DMFC,  which  takes  into  account  the  transport  phenomena  in 
all  the  layers  of  the  cell.  In  this  study,  a  two-dimensional  model 
of  a  FE-DMFC  has  been  developed  to  study  the  effects  of  the  fluid 
velocity  at  the  fuel,  air,  and  flowing  electrolyte  channel  inlets. 


transported  through  the  CBL.  The  electrochemical  reaction  occur¬ 
ring  at  the  CCL  is  shown  in  Eq.  (2). 

1 .502  +  6H+  +  6e"  -*  3H20  (2) 

The  methanol  crossing  over  the  membranes  and  FEC,  if  any,  also 
reaches  the  CCL.  Here,  the  methanol  reacts  with  oxygen  as  shown 
in  Eq.  (3). 

CH3OH  +  1 .502  -*  C02  +  2H20  (3) 

According  to  the  Ohm’s  law,  the  protonic  and  electronic  current 
densities  of  the  different  layers  can  be  shown  as  Eqs.  (4)  and  (5), 
respectively.  In  these  equations,  l  and  s  represent  the  electrolyte 
and  electrode  phases,  respectively. 

ix  =  -01 V0/  (for  ACL,  AM,  FEC,  CM,  and  CCL)  (4) 

is  =  — <ts V0s  (for  ABL,  ACL,  CCL,  and  CBL)  (5) 

The  charge  balance  for  the  electron  conduction  at  the  ABL  and 
CBL  is  shown  in  Eq.  (6);  and  the  charge  balance  for  the  proton 
conduction  at  the  AM,  FEC,  and  CM  is  shown  in  Eq.  (7).  Here,  it  is 
assumed  that  there  are  no  electron  or  proton  losses  in  these  layers. 


2.  Modeling 

A  two-dimensional  model  of  a  FE-DMFC  has  been  developed. 
For  validation  purpose,  we  have  also  developed  a  two-dimensional 
DMFC  model  using  similar  modeling  approach  and  equations  as 
the  FE-DMFC.  The  details  of  the  modeling  for  the  FE-DMFC  are  pre¬ 
sented  in  the  following  sub-sections.  The  main  assumptions  used 
in  the  modeling  are  as  follows. 

•  The  changes  of  the  protonic  and  electronic  conductivities  within 
a  layer  are  neglected. 

•  The  formations  of  the  C02  bubbles  and  water  vapor  are  neglected ; 
i.e.  two  phase  effects  are  not  taken  into  account. 

•  Membranes  are  fully  hydrated. 

•  Methanol  is  fully  consumed  at  the  interface  of  the  cathode  mem¬ 
brane  and  the  cathode  catalyst  layer. 

•  The  effects  due  to  the  channel  curvature  are  not  taken  into 
account. 

•  The  flow  in  the  electrolyte  channel  is  considered  as  fully  devel¬ 
oped  laminar  flow. 

•  The  fuel  cell  is  isothermal. 

•  The  fuel  cell  operates  at  the  steady  state  condition. 

2.2.  Proton  and  electron  transport  equations 

A  liquid  mixture  consisting  of  methanol  and  water  is  pumped  to 
the  fuel  cell  through  the  fuel  channel  (FC).  Some  amount  of  these 
species  diffuses  through  the  anode  backing  layer  (ABL),  which  is 
typically  carbon  cloth;  and  then  they  reach  the  anode  catalyst  layer 
(ACL)  having  Pt  and  Ru,  where  the  protons  and  electrons  are  gen¬ 
erated.  The  reaction  occurring  at  the  ACL  is  shown  in  Eq.  (1 ). 

CH3OH  +  H20  -*  C02  +  6H+  +  6e-  (1) 

The  protons  generated  in  the  ACL  reach  the  cathode  catalyst 
layer  (CCL)  having  Pt  after  they  transport  through  Nation®  anode 
and  cathode  membranes  (AM  and  CM)  and  flowing  electrolyte 
channel  (FEC)  consisting  of  diluted  sulfuric  acid  solution  (i.e.  H2S04 
and  H20).  The  electrons  are  conducted  to  the  external  load  through 
the  ABL  and  the  anode  interconnect.  These  electrons  reach  the  CCL 
after  they  pass  through  the  cathode  interconnect  and  the  cath¬ 
ode  backing  layer  (CBL).  Please  note  that  the  interconnects  are  not 
taken  into  account  in  this  study  due  to  their  high  electronic  con¬ 
ductivities.  As  the  electrons  and  protons  reach  the  CCL,  they  react 
with  oxygen  pumped  to  the  cell  through  the  air  channel  (AC)  and 


V  is  =  0 

(for  ABL  and  CBL) 

(6) 

V  ix  =  0 

(for  AM,  FEC,  and  CM) 

(7) 

The  charge  balances  for  the  proton  and  electron  conduction  at 
the  ACL  are  shown  in  Eqs.  (8)  and  (9),  respectively.  As  the  pro¬ 
tons  and  electrons  move  in  the  opposite  directions,  the  volumetric 
current  density  appears  as  a  source  term  for  the  proton  transport 
equation;  whereas  it  is  a  sink  term  for  the  electron  transport  equa¬ 
tion.  The  volumetric  current  density  of  the  anode  can  be  calculated 
using  Eq.  (10)  [20].  The  overpotential  at  the  anode  can  be  found 
using  Eq.  (11). 


V  •  it  =ja 

V  ■  is  =  ~ja 

•  =  aipf CMeQH  exp ((qaF/RT)r)a) 
QvieOH  +  Kc  exp ((aaF/RT)t]a) 

Va  =  <Ps  ~  <t>l  ~  EaQ 


(8) 

(9) 

(10) 

(11) 


The  charge  balances  for  the  proton  and  electron  conduction  at 
the  CCL  are  shown  in  Eqs.  (12)  and  (13),  respectively.  Due  to  the 
chemical  reaction  of  methanol  with  oxygen  (Eq.  (3))  occurring  in 
this  layer,  the  volumetric  current  density  of  the  cathode  decreases 
as  shown  in  Eq.  (14)  due  to  a  parasitic  loss  called  as  the  crossover 
current  density.  The  value  of  this  current  density  can  be  calculated 
using  the  average  molar  flow  rate  of  methanol  reaching  the  inter¬ 
face  of  CM  and  CCL  as  shown  in  Eq.  (15).  The  overpotential  at  the 
cathode  is  shown  in  Eq.  (16). 


V  if  =  -jc 

V  •  is  =  jc 


Jc  +Jxover  —  01 


ref  (-O2 


oc 


rref 

lo2 


6  FN 


Jxover  — 


//CM  |  CCL 
MeOH 


fcd 


exp 


Vc  =  4>s  -  <Pl  -  Ecq 


(12) 

(13) 

(14) 

(15) 

(16) 


The  boundary  conditions  applied  to  solve  the  electron  and  pro¬ 
ton  transport  equations  are  shown  in  Eqs.  ( 1 7)-(  1 9).  Eq.  ( 1 9)  shows 
the  insulation  boundary  condition  applied  to  the  remaining  bound¬ 
aries;  and  k  denotes  both  electrode  and  electrolyte  phases. 


0s,FC|ABL  =  0 


(17) 
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CH3OH+H2O+CO2  H2SO4+CH3OH+H2O  O2+N2+CO2+H2O 


0s,CBL|AC  =  Vcell  (18) 

-n  ik  =  0  (19) 

2.2.  Continuity  and  momentum  equations 

The  continuity  equation  for  FC,  ABL,  CBL,  and  AC  is  shown  in  Eq. 
(20);  and  this  equation  for  the  catalyst  layers  is  shown  in  Eq.  (21). 
In  the  latter  equation,  k  denotes  any  species,  e.g.  CH3OH,  H20,  02, 
and  C02,  produced  or  consumed  in  a  catalyst  layer.  For  this  layer, 
the  gain  or  loss  in  the  total  mass  of  the  fluid  mixture  originates 
from  the  production  or  consumption  of  the  protons.  The  source  or 
sink  terms  in  the  continuity  equation  of  the  catalyst  layers  can  be 
calculated  using  the  equations  derived  from  the  species  transport 
equations  as  discussed  in  Section  2.3. 


The  average  velocities  at  the  fuel  and  air  channel  inlets  are 
calculated  using  a  reference  condition  and  a  stoichiometric  flow 
coefficient  as  shown  in  Eqs.  (27)  and  (28),  respectively  [20].  In  these 
equations,  width  of  the  cell,  wce//,  is  equal  to  the  summation  of  the 
widths  of  the  channel  and  rib.  For  the  walls,  no-slip  condition,  Eq. 
(25),  is  applied. 


^abef^cell^cell 
vfc,in  ~  f 

®^MeOH  f/cW/c 

(27) 

^cbef^cell^ cell 

Vac, in  —  nr 

(28) 

4FCgCJntacwac 

2.3.  Species  transport  equations 

pV  u  =  0  (for  FC,  ABL,  CBL,  and  AC)  (20) 

pV  ■u  =  'V]sCikMWk  (for  ACL  and  CCL)  (21) 

The  momentum  equation  for  FC,  ABL,  ACL,  CCL,  CBL,  and  AC  can 
be  given  as  Eq.  (22)  [20].  For  the  AC  and  FC,  the  porosity  term  is 
equal  to  1.  In  this  equation,  the  source  term  is  zero  for  the  AC  and 
FC  as  shown  in  Eq.  (23);  whereas  its  value  can  be  calculated  using 
Darcy’s  law,  Eq.  (24),  for  the  ABL,  ACL,  CCL,  and  CBL.  For  the  AM 
and  CM,  the  fluid  velocity  is  assumed  to  be  zero  [20].  The  velocity 
distribution  in  the  FEC  is  found  considering  fully  developed  laminar 
flow  conditions  in  this  channel  as  shown  in  Eq.  (26)  [18]. 

4r(u  •  V)u  =  -  VP  +  Vt  +  Su  (for  FC,  ABL,  ACL,  CCL,  CBL,  and  AC) 

£z 

(22) 


2.3 A.  Methanol  transport 

Methanol  enters  the  fuel  cell  through  the  fuel  channel.  Due  to 
the  diffusion  and  convection  mechanisms,  some  portion  of  it  flows 
through  the  ABL  and  the  remaining  methanol  exits  through  the  fuel 
channel.  The  molar  flow  rate  of  methanol  and  the  transport  equa¬ 
tion  in  the  fuel  channel  are  shown  in  Eqs.  (29)  and  (30),  respectively. 


At"  =  -dvc  +  uc 

(29) 

V  •  (— DVC)  +  U  •  VC  =  0 

(30) 

In  the  ABL,  the  diffusion  and  convection  mechanisms  are  also 
effective  in  the  transport  of  methanol  to  the  ACL.  However,  an  effec¬ 
tive  diffusion  coefficient  should  be  used  considering  the  effects  of 
porosity  as  shown  in  Eqs.  (31 )— (33). 

N"  =  — DeffVC  +  uC  (31) 


Su  =  0  (for  AC  and  FC) 


Su  =  -^u  (for  ABL,  ACL,  CCL,  and  CBL) 


u  =  0  (for  AM  and  CM) 

6 Vfec(x  ~  *AM|FEC)(*FEC|CM  ~  x) 


Uy  — 


nrepeatWceutjec 


(for  FEC) 


(23)  V-  (-DeffVC)  +  u  VC  =  0  (32) 

(24)  where 

(25)  Deff  =  £1-5D  (33) 

In  the  ACL,  the  production  of  protons  leads  to  the  drag  of 

(26)  methanol  to  the  membrane  due  to  the  electro-osmosis  effect.  The 

molar  flow  rate  of  methanol  and  the  transport  equation  at  this  layer 
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can  be  shown  as  Eqs.  (34)  and  (35),  respectively.  Combining  Eq.  (35) 
and  (8),  the  transport  equation  can  be  shown  as  Eq.  (36). 

1 v"  =  -DeffVC  +  uC  +  n¥e0H  ^ 

a  p 


V  •  (-DeffVC)  +  u  •  VC  + 


V  •  (-DeffVC)  +  u  VC  = 


,MeOH 

ld 


In  Eqs.  (35)  and  (36),  the  drag  coefficient  of  diluted  methanol 
is  considered  as  constant  using  the  average  concentrations  of 
methanol  and  water  at  the  interface  of  ACL  and  AM,  which  is  shown 
in  Eq.  (37)  [20]. 


^MeOH  ~  „H20  OvieOH 

n“  =  ”*  -qT 


0 


(37) 


ACL|AM 


In  the  AM  and  CM,  the  convective  flow  is  neglected  as  mentioned 
in  Section  2.2.  Hence,  electro-osmosis  and  diffusion  are  the  only 
mechanisms  considered  as  shown  in  Eqs.  (38)  and  (39).  Combining 
Eqs.  (7)  and  (39),  the  transport  equation  reduces  to  Eq.  (40). 


crossover  and  the  electrochemical  reaction  appear  as  a  sink  term. 
The  molar  flow  rate  and  transport  equation  for  oxygen  at  the  CCL 
are  shown  in  Eqs.  (47)  and  (48),  respectively. 


(34) 

'  // 

N  = 

V  •  h  -ja 

~  6  F 

(35) 

V  •  (  — 

'  i  „MeOH\ 

6+n<>  ) 

(36) 

2.3.4. 

V  •  (—  DeffVC)  +  U  •  VC  —  — —  (jc  -fjxover) 

4  r 


(47) 


(48) 


The  concentration  of  the  methanol  at  the  fuel  channel  inlet  is 
known,  Eq.  (49);  and  the  concentrations  of  the  water  at  the  fuel 
channel  inlet  and  oxygen  at  the  air  channel  inlet  can  be  calculated 
from  the  given  input  data.  As  the  fuel  storage  tank  of  the  fuel  cell 
system  consists  of  water  and  methanol  with  a  constant  volume, 
the  concentration  of  water  at  the  fuel  channel  inlet  can  be  given 
in  terms  of  the  concentration  of  methanol  as  shown  in  Eq.  (50). 
The  oxygen  concentration  at  the  inlet  of  the  air  channel  can  be 
expressed  in  terms  of  the  molar  ratio  of  oxygen  in  the  air  and  inlet 
pressure  of  the  air  channel  as  shown  in  Eq.  (51). 


^MeOH  ~~  Known 


(49) 


n"  =  -DeffVC +  n“e0H^ 

nMeOHv  .  i 

V  •  (-DeffVC)  H — - — — - -  =  0 

V- (-DeffVC)  =  0 


(38) 

(39) 

(40) 


ffc,in  _  Ph2Q 
h2°  MWh2o 


CMeOHMVVMeOH 


PMeOH 


P 

rac,m  _  v  1  ac,in 

lo2  -  *o2  RT 


(50) 

(51) 


In  the  FEC,  the  effect  of  the  convective  flow  of  sulfuric  acid  solu¬ 
tion  should  be  taken  into  account.  Using  a  similar  approach  with  a 
previous  study  conducted  by  Colpan  et  al.  [18],  this  effect  has  been 
considered  as  a  sink  term  in  the  methanol  transport  equation.  The 
molar  flow  rate  of  methanol  is  shown  in  Eq.  (41 ).  Using  the  velocity 
of  the  FEC  shown  in  Eq.  (26),  the  methanol  transport  equation  at 
this  layer  can  be  shown  as  Eq.  (42).  Combining  Eqs.  (7)  and  (42),  the 
methanol  transport  equation  reduces  to  Eq.  (43). 

n"  =  -DeffVC  +  n“e0H  lj  (41 ) 

V  (_ d  ffVC)  +  l1d1e°HV  1  h  =  S^/ecC^-XAMIFEcXXFEQCM  ~X)(12) 

F  II  repeat  C:eIlwceH  Jec 


In  this  study,  it  is  assumed  that  the  entire  methanol  reaching 
the  CCL  is  consumed  at  the  interface  of  CM  and  CCL  as  shown  in  Eq. 
(52).  It  should  be  noted  that  if  the  flow  rate  of  sulfuric  acid  is  high 
enough,  the  concentration  of  methanol  can  become  zero  within  the 
FEC  or  CM. 


rCM|CCL  _  n 
LMeOH  ~~  U 


(52) 


Continuity  equations  are  used  for  all  the  interior  boundaries.  The 
outflow  boundary  condition,  Eq.  (53),  for  the  channel  outlets  and 
no-flux  boundary  condition,  Eq.  (54),  for  the  remaining  boundaries 
are  applied. 


-it  •  VCfc  =  0 


(53) 


V-  (-DeffVC)  =  - 


6VfecC(x  -  Xam|Fec)(xfeqcm  -  x) 

3 

nrepeatLceiiWcentjrec 


(43) 


2.3.2.  Water  transport 

The  main  equations  for  water  transport  at  the  FC  and  ABL  are 
similar  with  those  for  methanol  transport  given  in  Eqs.  (29)-(33). 
For  the  ACL,  diffusion,  convection,  and  electro-osmosis  affect  the 
transport  of  water  as  shown  in  Eqs.  (44)  and  (45).  Combining  Eqs. 
(8)  and  (45),  the  transport  equation  can  be  shown  as  Eq.  (46).  In 
this  study,  we  assumed  a  constant  concentration  profile  for  water 
at  the  AM,  FEC,  and  CM  [18]. 

n"  =  -DeffVC  +  uC  +  n^°j 

nf°V  ■  i,  -j 
V-  (-DeffVC)  +  u-VC+-d  Ja 


F 

~ ja  (  1 


6F 


V  ■  (-DeffVC)  +  u  .  VC  =  -p  (g  +  nf° 


2.3.3.  Oxygen  transport 

The  main  equations  for  the  oxygen  transport  at  the  AC  and 
CBL  are  similar  to  those  for  the  methanol  transport  equations,  Eqs. 
(29)-(33).  In  the  CCL,  the  oxygen  consumption  due  to  the  methanol 


n  Nk  =  0 


2.4.  Output  parameters 


(54) 


The  average  current  density  of  the  cell  can  be  calculated  using 
Eq.  (55). 


^cell  ~ 


jadxdy 


(55) 


ACL 


The  power  density  of  the  cell,  and  the  electrical  efficiency  of  the 
sing  Eqs.  (56)  and  (57),  respectively. 

(56) 

(57) 


(44) 

K„ 

—  icell  ^ cell 

(45) 

4el  = 

(46) 

NMeOHLHV 

If  the  flowing  electrolyte  is  recirculated,  the  electrical  efficiency 
can  be  shown  as  follows. 


4el  = 


Wce,l 


(N, 


ii"fc,in  ^  fee,  out 


MeOH  MeOH 


)LHV 


(58) 
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Fig.  3.  Validation  of  the  2D  DMFC  model. 

3.  Results  and  discussion 

The  results  and  discussion  for  the  validation  of  the  model  and 
the  parametric  studies  including  the  effects  of  the  fluid  velocities 
at  fuel,  air,  and  flowing  electrolyte  channel  inlets  are  presented  in 
this  section.  The  main  input  data  used  in  the  simulations  are  shown 
in  Table  1.  The  solution  of  the  set  of  the  equations  given  in  Section 
2  is  done  using  Comsol  Multiphysics  4.2,  which  is  software  based 
on  the  finite  element  method.  After  a  grid  independent  analysis, 
51,817  triangular  elements  are  taken  for  the  geometry  chosen. 

3.2.  Validation  of  the  model 

The  results  of  the  2D  DMFC  model  developed  are  compared  to 
the  experimental  data  published  by  Ge  and  Liu  [21  ].  The  geometri¬ 
cal  dimensions  of  the  cell  and  the  main  operating  parameters  are 
given  in  Table  1 .  In  their  experiments,  the  following  materials  were 
used  as  the  components  of  the  cell:  Nation®  1 1 7  for  the  membrane, 
carbon  cloth  for  the  anode  and  cathode  backing  layers,  Pt-Ru  with 
a  loading  3  mg  cm-2  for  the  anode  catalyst  layer,  and  Pt-black  with 
a  loading  of  3  mg  cm-2  for  the  cathode  catalyst  layer.  The  compari¬ 
son  of  the  experimental  data  and  the  results  of  the  simulation  of  the 
model  between  the  cell  voltages  of  0.2  V  and  0.8  V  is  shown  in  Fig.  3. 
As  can  be  seen  from  this  figure,  there  is  a  good  agreement  between 
the  experimental  and  modeling  data.  The  minor  discrepancy  at  the 
low  current  density  conditions  can  be  attributed  to  the  simplified 
equations  used  for  electrochemical  relations,  Eqs.  (10)  and  (14), 
occurring  at  the  anode  and  cathode  catalyst  layers.  At  high  current 
densities  such  as  3000  A  m-2,  the  experimental  results  show  a  devi¬ 
ation  from  linearity  due  to  the  diffusion  limitations ;  however  as  this 
deviation  occurs  at  higher  current  densities  (i.e.  current  densities 
higher  than  3000 Am-2)  for  the  model  developed,  the  modeling 
results  show  a  linear  trend  for  the  high  current  density  conditions 
shown  in  Fig.  3.  The  reason  of  this  discrepancy  at  these  conditions 
can  be  attributed  to  neglecting  the  two-phase  effects  in  the  model 
developed;  which  has  a  significant  effect  on  the  limiting  current 
density. 

3.2.  Effect  of  the  fluid  velocity  at  the  fuel  channel  inlet 

The  velocity  of  the  diluted  methanol  solution  at  the  fuel  channel 
inlet  can  be  represented  by  the  stoichiometric  flow  coefficient  as 
shown  in  Eq.  (27).  For  the  given  coefficient,  methanol  concentra¬ 
tion  at  the  inlet,  geometry  of  the  cell,  and  chosen  reference  current 
density,  which  is  taken  as  1500 Am-2,  this  velocity  is  first  calcu¬ 
lated.  In  the  simulations,  the  volumetric  flow  rate  of  the  sulfuric 


Table  1 

Input  data  for  the  base  case. 


Input  parameter 

Value 

Length  of  the  cell 

6.5  x  10“2  m  [8] 

Width  of  the  air  and  fuel  channels 

1  x  10-3  m  [8] 

Width  of  the  repeat  element 

2  x  10“3m  [8] 

Thickness  of  the  air  and  fuel 

8  x  10-4  m  [8] 

channels 

Thickness  of  the  anode  and 

1.4  x  10“4  m  [8] 

cathode  backing  layers 

Thickness  of  the  anode  and 

3  x  10“5m  [8] 

cathode  catalyst  layers 

Thickness  of  the  membrane 

1.83  x  10“4m  [22] 

(DMFC) 

Thickness  of  the  anode  and 

9.2  x  10“5m  [18] 

cathode  membranes  (FE-DMFC) 
Thickness  of  the  flowing 

5  x  10“4m  [18] 

electrolyte  channel 

Number  of  repeat  elements  in  a 

32 

single  cell 

Molar  concentration  of  the 

2000  mol  lrr3 

methanol  at  the  fuel  channel 
inlet 

Molar  ratio  of  oxygen  at  the  air 

21% 

channel  inlet 

Stoichiometric  flow  coefficient  at 

2 

fuel  channel  inlet 

Stoichiometric  flow  coefficient  at 

3 

air  channel  inlet 

Volumetric  flow  rate  of  the  flowing 

1.67  x  10“7  m3  s”1 

electrolyte  per  single  cell 
Temperature  of  the  cell 

70  °C 

Pressure  of  the  anode  and  cathode 

1  atm 

outlets 

Porosity  of  the  anode  and  cathode 

0.6  [8] 

backing  layers 

Porosity  of  the  anode  and  cathode 

0.4  [8] 

catalyst  layers 

Porosity  of  the  membranes 

0.28  [8] 

Porosity  of  the  spacer  in  the 

0.6  [18] 

flowing  electrolyte  channel 
Permeability  of  the  anode  and 

2  x  10“12m2  [23] 

cathode  backing  layers 

Permeability  of  the  anode  and 

10“13m2  [23] 

cathode  catalyst  layers 

Permeability  of  the  flowing 

2  x  10“12  m2  (assumed) 

electrolyte  channel 

Coefficient  of  diffusion  of  methanol 

2.8  X  10-9  X  e[2436((l/353)-(l/T))]  m2  g-l 

in  water 

[24] 

Coefficient  of  diffusion  of  methanol 

4.9  X  10-10  X  e[2436((l/333)— (l/T))]  m2  S_1 

in  Nafion® 

[24] 

Coefficient  of  diffusion  of  oxygen  in 

(T  27J72X™  )m2s-'[25] 

gas 

Electro-osmotic  drag  coefficient  of 

1.6767  +  0.0155  x  T+  8.9074  x  10“5  x  T2 

water 

(T  is  in  °C)  [26] 

Cell  voltage 

0.4  V 

Protonic  conductivity  of  the 

10  S  m-1  [22] 

membrane 

Protonic  conductivity  of  the 

1455m-1  [18,27] 

sulfuric  acid  solution 

Electronic  conductivity  of  the 

300  5  m-1  [28] 

backing  layers 

Methanol  oxidation  reaction 

2.265  x  10-3  mol  m-3  [29] 

constant 

Reference  oxygen  concentration 

0.472  molm-3  [8] 

Anodic  transfer  coefficient 

0.5  [8] 

Cathodic  transfer  coefficient 

0.5  [8] 

Reference  exchange  current 

6  x  105  Am-3  (calibrated) 

density  times  specific  area  at 
anode 

Reference  exchange  current 

200 Am-3  [8] 

density  times  specific  area  at 
cathode 

Reference  current  density  for 

1500 Am-2  [23] 

calculating  the  fluid  velocity 
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Fig.  4.  Methanol  concentration  distribution  of  the  FE-DMFC  for  different  fuel  channel  inlet  velocities:  (a)  |a  = 1  or  VfCjn  =2.11x1 0-4  m  s-1 ,  (b)  |a  =  2  or  VfCjn  =  4.21  x  1 0-4  m  s_1 , 
(c)  |a  =  3  or  VfCjn  =  6.32  x  1 0-4  m  s_1 ,  and  (d)  |a  =  4  or  i ?fcin  =  8.42  x  1 0-4  m  s_1 . 


acid  solution  is  taken  high  enough;  which  means  that  the  amount 
of  methanol  crossover  to  the  cathode  side  is  negligible  for  all  cases. 

The  effect  of  the  fuel  channel  inlet  velocity  on  the  methanol 
concentration  distribution  for  different  stoichiometric  flow  coef¬ 
ficients,  at  the  cell  voltage  of  0.4  V,  is  shown  in  Fig.  4a-d.  It  can 
be  seen  from  these  figures  that,  for  any  given  flow  rate,  methanol 
concentration  decreases  along  the  thickness  and  flow  directions. 
The  main  reasons  for  these  decreases  are  the  consumption  of 
methanol  at  the  ACL  due  to  the  electrochemical  reaction  and  the 
drag  of  methanol  through  AM  due  to  the  electro-osmosis  effect.  The 
crossover  methanol  is  almost  completely  washed  away  at  the  FEC 
due  to  the  convective  effect  of  the  sulfuric  acid  solution  pumped 
to  the  cell.  The  remaining  little  amount  of  methanol  is  spent  dur¬ 
ing  the  chemical  reaction  occurring  at  the  CCL  In  this  study,  we 
assume  that  this  reaction  occurs  at  the  interface  CCL  and  CM; 
and  the  concentration  of  methanol  becomes  zero  at  this  interface. 
Fig.  4a-d  also  shows  that  as  the  stoichiometric  flow  coefficient  or 
the  fuel  channel  inlet  velocity  increases,  the  concentration  differ¬ 
ence  between  the  inlet  and  exit  of  the  cell  decreases.  For  example, 
the  average  methanol  concentration  at  the  exit  of  the  fuel  channel 
is  828  mol  m-3,  1353  mol  nrr3,  1547  mol  nrr3,  and  1649  mol  nrr3 


when  the  stoichiometric  flow  coefficient  at  the  fuel  channel  inlet  is 
taken  as  1,  2,  3  and  4,  respectively,  for  the  given  input  data.  These 
results  show  that  more  methanol  is  wasted  (i.e.  less  methanol  is 
utilized)  when  the  fluid  velocity  at  the  fuel  channel  inlet  increases. 

The  effect  of  the  fluid  velocity  at  the  inlet  of  the  fuel  channel 
on  the  cell  voltage  and  power  density  is  shown  in  Fig.  5.  As  it  can 
be  seen  from  this  figure,  the  trend  of  the  polarization  curves  does 
not  change  significantly  for  different  stoichiometric  coefficients. 
However,  the  limiting  current  density  increases  as  this  coefficient 
increase;  thus  lower  fuel  channel  inlet  velocities  cause  a  more  lim¬ 
ited  operating  range  for  the  fuel  cell.  These  trends  are  compared 
with  the  experimental  study  by  Ge  and  Liu  [21  ];  and  it  is  seen  that 
they  are  in  good  agreement  with  each  other.  In  addition,  this  fig¬ 
ure  shows  that  higher  maximum  power  densities  could  be  achieved 
with  higher  stoichiometric  coefficients.  Numerically,  the  maximum 
power  densities  are  found  as  600  Wm~2,  670 Wm-2,  720  Wm-2, 
and  730 Wm-2  when  the  stoichiometric  flow  coefficients  at  the 
fuel  channel  inlet  are  taken  as  1,  2,  3  and  4,  respectively.  On  the 
other  hand,  an  increase  in  the  flow  rate  of  the  methanol  pumped 
to  the  fuel  cell  causes  an  increase  in  the  power  demand  of  the  fuel 
pump;  however  this  demand  is  not  considered  in  this  study  as  the 
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Fig.  5.  Cell  voltage  and  power  density  of  the  FE-DMFC  for  different  fuel  channel  inlet 
velocities.  (The  limiting  current  densities  are  shown  with  arrows.) 
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Fig.  6.  Electrical  efficiency  of  the  FE-DMFC  and  recirculated  FE-DMFC  for  different 
fuel  channel  inlet  velocities. 
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Fig.  7.  Oxygen  concentration  distribution  of  the  FE-DMFC  for  different  air  channel  inlet  velocities:  (a)  £c  =  1  or  vac,m  =  8-47  x  1 0-2  m  s-1 ,  (b)  £c  =  2  or  vac,m  =  1 6-94  x  1 0-2  m  s-1 
(c)  I'c  =  3  or  vac, in  =  25.41  x  1 0-2  m  s-1 ,  and  (d)  £c  =  4  or  vac,m  =  33.88  x  1 0-2  m  s_1 . 
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focus  of  the  study  is  the  cell  itself  rather  than  the  complete  system 
including  balance  of  plant. 

The  fluid  velocity  at  the  fuel  channel  inlet  affects  the  electrical 
efficiency  of  the  cell  significantly,  as  can  be  seen  in  Fig.  6.  The  ratio 
of  the  total  amount  of  the  methanol  utilized  to  the  methanol  enter¬ 
ing  the  fuel  cell  increases  as  the  fluid  velocity  decreases;  which 
in  turn  leads  to  an  increase  in  the  electrical  efficiency  of  the  cell. 
Hence,  lower  stoichiometric  coefficients  at  the  fuel  channel  inlet 
are  desirable  to  obtain  higher  electrical  efficiencies.  However,  it 
should  be  kept  in  mind  that  lower  maximum  power  densities  are 
also  achieved  at  these  coefficients.  Therefore,  since  power  density 
and  electrical  efficiency  cannot  be  maximized  at  the  same  time,  the 
choice  of  the  fluid  velocity  at  the  fuel  channel  inlet  mainly  depends 
on  the  purpose  of  the  device  and  economical  factors  such  as  the  fuel 
cost  and  the  purchase  and  operation  costs  of  the  fuel  cell.  However, 
these  considerations  are  out  of  the  scope  of  the  paper. 

The  electrical  efficiency  of  the  FE-DMFC  can  be  further  increased 
if  the  methanol  leaving  the  FEC  can  be  recovered,  separated  from 
sulfuric  acid  and  recirculated  back  to  the  methanol  storage  tank. 
Kordesch  et  al.  [7]  proposed  such  a  configuration  including  a  sep¬ 
arator  to  their  basic  design.  As  less  amount  of  methanol  enters  the 
fuel  cell  at  a  given  time  due  to  the  addition  of  this  recirculation 
loop,  the  electrical  efficiency  of  a  recirculated  FE-DMFC  at  a  given 
stoichiometric  coefficient  becomes  higher  than  that  of  a  FE-DMFC, 
which  can  also  be  seen  in  Fig.  6.  The  maximum  electrical  efficiency 
that  can  be  achieved  for  the  recirculated  FE-DMFC  is  found  to  be  56% 
when  the  stoichiometric  coefficient  is  taken  as  1  and  the  current 
density  is  equal  to  1700  Am-2. 

3.3.  Effect  of  the  fluid  velocity  at  the  air  channel  inlet 

The  velocity  of  the  air  entering  the  fuel  cell  is  another  impor¬ 
tant  parameter  that  can  affect  the  performance  of  the  system.  This 
velocity  can  be  represented  by  a  stoichiometric  flow  coefficient  as 
shown  in  Eq.  (28).  The  value  of  this  velocity  depends  on  this  coef¬ 
ficient,  the  geometry  of  the  cell,  oxygen  concentration  at  the  air 
channel  inlet,  and  the  reference  current  density,  which  is  taken  as 
1 500  A  m- 2 .  In  this  paper,  the  effect  of  this  coefficient  on  the  oxygen 
concentration,  cell  voltage,  power  density,  and  electrical  efficiency 
is  studied.  On  the  other  hand,  this  coefficient  has  also  an  effect  in 
preventing  the  cathode  flooding  and  the  power  demand  of  the  air 
blower;  however  these  effects  are  not  considered  in  this  paper. 

The  oxygen  consumption  at  the  CCL  due  to  the  electrochemical 
reaction  and  the  oxidation  of  the  crossover  methanol  results  in  a 
decrease  in  the  oxygen  concentration  at  the  flow  direction,  which 
can  be  seen  in  Fig.  7a-d.  These  figures  also  show  the  oxygen  con¬ 
centration  distribution  in  a  FE-DMFC  for  different  stoichiometric 
flow  coefficients  when  the  cell  voltage  is  0.4  V.  The  concentration 
difference  along  the  inlet  and  exit  decreases  as  this  coefficient 
increases.  Numerically,  for  the  oxygen  concentration  at  the  air 
channel  inlet  as  7.46  mol  m-3,  this  concentration  at  the  air  channel 
outlet  is  calculated  as  4.45  mol  m-3,  5.85  mol  m~3,  6.35  mol  m-3, 
and  6.61  mol  m~3  when  the  stoichiometric  flow  coefficients  at  the 
air  channel  inlet  are  taken  as  1,2,  3,  and  4,  respectively.  These 
results  show  that  more  air  is  wasted  (i.e.  less  air  is  utilized)  for 
higher  coefficients. 

The  effect  of  the  air  velocity  entering  the  cell  on  the  cell  voltage 
and  power  density  of  the  FE-DMFC  is  shown  in  Fig.  8.  This  figure 
shows  that  there  is  a  slight  increase  in  the  performance  of  the  cell 
as  the  velocity  increases.  The  trend  of  this  increase  is  in  good  agree¬ 
ment  with  the  experimental  data  published  by  Ge  and  Liu  [21  ].  This 
increase  can  be  explained  as  follows:  the  cathodic  current  density 
is  proportional  to  the  concentration  of  oxygen  at  the  CCL  as  shown 
in  Eq.  (14).  This  concentration  increases  with  an  increase  in  the 
air  velocity;  hence  the  cathodic  overpotential  decreases  and  the 
performance  of  the  cell  increases. 


Cun  eut  Density  [A  in  2l 

Fig.  8.  Cell  voltage  and  power  density  of  the  FE-DMFC  for  different  air  channel  inlet 
velocities. 


Fig.  9  shows  the  electrical  efficiency  of  the  cell  for  the  FE-DMFC 
and  recirculated  FE-DMFC  for  different  stoichiometric  flow  coeffi¬ 
cients  at  the  air  channel  inlet.  This  figure  shows  that  as  the  air  flow 
rate  increases  the  electrical  efficiency  increases  since  there  is  no 
change  in  the  fuel  entering  the  system  and  power  density  increases. 
As  discussed  in  Section  3.2,  this  efficiency  can  be  further  increased 
if  a  recirculation  loop  is  used  for  the  methanol  leaving  the  FEC.  For 
example,  the  electrical  efficiencies  of  the  recirculated  FE-DMFC,  for 
the  cell  voltage  as  0.4  V  and  stoichiometric  flow  coefficient  at  the 
fuel  channel  inlet  as  2,  are  found  as  17.8%,  18.8%,  19.1%,  and  19.2% 
when  the  stoichiometric  flow  coefficients  at  the  air  channel  inlet 
are  taken  as  1 , 2, 3,  and  4,  respectively.  Here,  it  should  be  noted  that 
the  trend  of  the  electrical  efficiency  of  the  system  including  balance 
of  plant  such  as  air  blower  could  be  slightly  different  as  the  power 
demand  of  the  blower  increases  slightly.  However,  this  demand  is 
not  considered  in  this  study. 

3.4.  Effect  of  the  fluid  velocity  at  the  flowing  electrolyte  channel 
inlet 

The  amount  of  methanol  reaching  the  cathode  side  can  be 
mainly  controlled  by  adjusting  the  velocity  of  the  sulfuric  acid 
pumped  through  the  FEC.  Assuming  a  fully  developed  laminar  flow 
condition,  this  velocity  varies  only  with  thickness  direction  and 
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Fig.  9.  Electrical  efficiency  of  the  FE-DMFC  and  recirculated  FE-DMFC  for  different 
air  channel  inlet  velocities. 
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Fig.  10.  Crossover  current  density  of  the  FE-DMFC  for  different  volumetric  flow  rate 
of  flowing  electrolyte. 

depends  on  the  volumetric  flow  rate  of  the  sulfuric  acid,  num¬ 
ber  of  repeat  element,  and  the  geometry  of  the  cell,  as  shown  in 
Eq.  (26).  Fig.  10  shows  the  change  of  the  volumetric  crossover 
current  density,  which  is  a  measure  of  the  amount  of  crossover 
methanol,  with  different  volumetric  flow  rates  of  sulfuric  acid: 
0.1  ml  min-1  (1.67  x  10-9  m3  s_1),  1  ml  min-1  (1.67  x  10~8  m3  s_1), 
and  10  ml  min-1  (1.67  x  10~7  m3  s_1).  These  flow  rates  corre¬ 
spond  to  the  following  average  fluid  velocities  respectively: 
5.21  x  10-5  ms-1,  5.21  x  lO^ms"1,  and  5.21  x  lO^ms”1.  This 
figure  shows  that  high  flow  rates  such  as  10  ml  min-1  should  be 
taken  to  eliminate  the  effect  of  the  methanol  crossover.  This  result 
is  as  expected  since  the  convective  flow  at  the  flow  direction  gets 
more  influential  to  sweep  away  the  crossover  methanol  at  the 
higher  flow  rates. 

The  effect  of  the  volumetric  flow  rate  of  the  sulfuric  acid  on 
the  cell  voltage  and  power  density  of  the  FE-DMFC  is  shown  in 
Fig.  11.  In  a  FE-DMFC,  the  addition  of  the  FEC  increases  the  ohmic 
losses;  whereas  it  decreases  the  cathodic  overpotential.  When  the 
thickness  of  the  FEC  is  kept  at  its  minimum  value,  the  ohmic  losses 
associated  with  this  FEC  could  be  minimized;  and  the  flow  rate  of 
the  FEC  becomes  the  main  operating  parameter  affecting  the  per¬ 
formance  of  the  cell.  Figs.  11  and  12  show  that  the  performance 
of  the  cell  slightly  increases  with  an  increase  in  the  flow  rate  of 
the  FEC.  However,  the  electrical  efficiency  of  the  cell  increases 
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Fig.  11.  Cell  voltage  and  power  density  of  the  FE-DMFC  for  different  volumetric  flow 
rate  of  flowing  electrolyte. 
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Fig.  12.  Electrical  efficiency  of  the  FE-DMFC  and  recirculated  FE-DMFC  for  different 
volumetric  flow  rate  of  flowing  electrolyte. 


substantially  if  a  recirculation  loop  is  used,  as  shown  in  Fig.  12. 
For  example,  when  the  stoichiometric  flow  coefficients  at  the  fuel 
and  air  channel  inlets  are  taken  as  2  and  3,  the  maximum  electrical 
efficiency  and  power  density  that  can  be  achieved  are  found  as  26% 
and  670  Wm-2,  respectively. 


4.  Conclusions 

A  two-dimensional  multi-physics  model  has  been  developed  to 
simulate  the  performance  of  the  cell  and  to  study  the  effects  of  some 
of  the  key  operating  parameters.  These  parameters  include  the  fluid 
velocity  at  the  fuel,  air,  and  flowing  electrolyte  channel  inlets.  Using 
this  model,  the  concentration  distributions  of  the  species,  such  as 
methanol  and  oxygen,  along  the  thickness  and  flow  directions,  the 
cell  voltage,  the  power  density,  and  the  electrical  efficiency  of  the 
cell  are  found.  The  main  findings  of  this  study  are  listed  below. 

•  There  is  a  good  agreement  between  the  results  of  the  model  and 
the  experimental  data. 

•  As  the  stoichiometric  flow  coefficient  at  the  fuel  channel  inlet 
increases  (i.e.  as  the  fluid  velocity  at  the  fuel  channel  inlet 
increases),  higher  power  densities  could  be  achieved;  however 
less  methanol  is  utilized  and  the  electrical  efficiency  of  the  cell 
decreases. 

•  An  increase  in  the  stoichiometric  flow  coefficient  at  the  air  chan¬ 
nel  (i.e.  an  increase  in  the  fluid  velocity  at  the  air  channel  inlet) 
yields  to  an  increase  in  both  the  power  density  and  the  electrical 
efficiency  of  the  cell. 

•  The  electrical  efficiency  of  the  FE-DMFC  can  be  further  increased 
if  the  methanol  leaving  the  FEC  is  recirculated  to  the  methanol 
storage  tank. 

•  The  velocity  at  the  FEC  inlet  should  be  taken  high  enough  to 
eliminate  the  methanol  crossover  to  the  cathode  side. 

In  this  study,  some  operating  parameters  are  suggested  to 
increase  the  performance  of  a  FE-DMFC.  The  results  and  main 
findings  of  the  simulations  conducted  are  expected  to  help  in  the 
development  of  the  FE-DMFC  to  be  used  in  portable  applications.  A 
3D  model  will  be  developed  as  a  future  work  to  have  a  more  accu¬ 
rate  model  and  to  find  the  distributions  of  the  output  parameters 
in  all  three  directions. 
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